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Abstract: Fully coupled simulations of hydrodynamics and radiative transfer are essential to a number of fields 
ranging from astrophysics to engineering applications. Of particular interest in this work are hypersonic 
atmospheric entries and associated experimental apparatus, e.g., shock tubes and high enthalpy testing facilities. 
The radiative transfer calculations must supply to the CFD a heating term in the energy equation in the form of 
the divergence of the radiative heat flux and the radiative heat fluxes to bounding surfaces. It is most efficient to 
solve the radiative transfer equation on the same grid as the CFD solution, and this work presents an algorithm 
with improved accuracy for such simulations on structured and unstructured grids compared to more 
conventional approaches. Results will be shown for shock radiation during hypersonic reentry. Issues of 
parallelization within a radiation sweep will also be discussed. 
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1 Introduction 

The radiative transfer equation (RTE) for the radiative intensity I in direction Q in a non- scattering medium is given by 


Q-V/(x,Q)-*(x)(S(x)-/(x,Q)) 


( 1 ) 


where S is the source function and x is the opacity. Frequency dependence has been suppressed for simplicity of notation. 
Upon integration over a computational cell and discretizing, approximating the cell volume integral with the values at the 
center, one obtains 

2 I(x t ,Q)Q-AS t -x(\)(S(\)-I(x c ,Q))V (2) 

faces k 


where x c is the position of the cell center, x k is the center of face k , V is the cell volume, and AS*. is the surface area vector 
of face k of the cell. Splitting the surface summation into incoming and outgoing parts and using cell-center values for the 
outgoing facial intensities in each cell, one gets the standard finite- volume method for radiative transfer [1]: 


I(x c ,Q) = 


2 /(x„Q)|Q-AS,| + x(x c )S(x c )V 

&,Q-AS*<0 

2 Q - AS * + *( x c)r 

k,Q-AS k >0 


( 3 ) 


2 Improved formulation 

The relatively inaccurate volume-integral approximation in (2) and the approximation of using cell-center values of I at the 
outgoing faces in (3) can both be eliminated by using the 1-d exact radiative transfer solution for constant x and S along 
direction Q. To determine the effective cross-cell propagation distance d in direction Q , we divide the volume V by the 
total area of the incoming (or, equivalently, outgoing) faces projected on that direction: 

d ~ 2 |Q-AS,| <4) 

i\Q- A$<0 


Then the intensity on each outgoing face o is given by 



( 5 ) 


2 /(x,.,Q)|Q-AS,| 

/(x„,Q) = l ' n " S ' < V Ir, a C I e ~ Xt '' >d + (1 - <’‘* (X,W >' S W 

2j • AS f | 

/,Q-ASj<0 

where x Q indicates an outgoing face center and x, an incoming face center. This method accounts more accurately for 
emission and absorption acting between the incoming and outgoing faces of a computational cell, essentially by treating 
these processes with sub-cell accuracy using an analytic form. Eqs. (5) and (3) are the same to first order in the optical 
thickness of the cell, x( x c)^« Values of I are only required at the cell faces. The total radiative heating in a cell, needed for 
the CFD energy equation, is easily computed as 

Ja 2 (6) 

faces k 

No differences are needed to compute this quantity, only a summation over the cell faces. The other quantity typically 
needed in a rad-hydro calculation is the heat flux at a body surface. Assuming that face k of a cell lies on this surface, the 
heat flux through that surface is simply: 

q* ** f a toi(*k> Q ) dQ ( ? ) 


3 Verification and Comparison 

Figure 1 below shows a comparison between a Cartesian-mesh, short-characteristic-based RTE solution (leftmost), the 
same problem on a pyramidal mesh using the method described here (middle), and a solution using the classical method 
(eq. (3)) on the same pyramidal mesh (rightmost). Contours of cell-average div(q) are shown in the midplane for a 
hyperboloidal radiative source emulating a reentry shock around a bluff body; x is streamwise and y is cross-stream. The 
rightmost plot shows a loss of heating in the classical method compared to the other two, presumably due to 
approximations in eq. (3) that have been improved using the current method (eq. (5)). The characteristic method result 
shows it to be slightly more diffusive, as expected, compared to both finite- volume methods, but the approximate 
agreement serves to verify the implementations of the latter. The final paper will make extensive comparisons for media of 
various degrees of optical thickness to show the advantages of the sub-cell treatment of emission and absorption. 



Fig.l: Comparison of div(q) for Cartesian short-characteristic, and current and classical finite volume methods. 
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